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Abstract 

The purpose of this note is to give the one loop radiative corrections to the 
top quark pair production in the pp annihilation at the Fermilab Tevatron 
in the context of the Minimal Supersymmetric Model. We concentrate here 
on the supersymmetric QCD corrections and give the analytic expression for 
these corrections. Recently Li et. al. have reported the supersymmetric QCD 
corrections to this process we indicate clearly a comparison of their and our 
work. In particular, we find additional corrections [crossed box and gluon self- 
energy] at the one loop level which are not given by Li et. al.. Our numerical 
results disagree with the original claim of Li et. al. The numerical values given 
by them in a recent erratum do agree with the general trend of our numerical 
results however the actual values still disagree. We find that the percentage 
corrections at the hadronic corrections changes from 22% to —0.5% as the 
squark mass is changed from 100 GeV to 600 GeV, for a gluino mass of 200 
GeV. For a gluino mass of 150 GeV the squark-mass dependence is less abrupt; 
they change from —5.3% to 1% as the squark mass is varied between 100 GeV 
and 600 GeV. We also present numerical results for differential cross section 
at the hadronic level, and percentage corrections at the parton level. 
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I. INTRODUCTION 



As is well known by now the top quark existence has been experimentally shown by 
the CDF |T[]and DO [0 at almost 100% confidence level. Two interesting parameters, the 
mass of top and the cross section for top pair production, have been found as follows: by 
the CDF i 

1. mf"*- = 176 ± 9 GeV, 

The DO finds for the same parameters 

1. mfP*- = 170 ± 18 GeV, 

2. cT^'f*- =5.2 ±1.8 pb. 

The standard model theoretical predictions for the top pair production cross section is, 
assuming a top mass of 170 and 175 GeV 0], 

1. al^'"'' = 6.481°:°^ pb, mt = 170 GeV, 

2. cr*f = 5.521°:°^ pb, = 175 GeV. 

A theoretical fit based on the Standard Model Electroweak Precision calculations gives for 
the top mass the following limits 0] 

mt = 179 ± rtf^tlZTcS^r T 2(«.) T 5(5J (1) 

here as = 0.120 ± 0.07 and 5a = 0.03 ± 0.09 are the uncertainities in as(mz) and tt(m|) 
M, respectively. 



Once the main injector upgrade becomes operational in 1999 ^ at Fermilab, the exper- 
imental sensitivity will be highly increased. For example the uncertainty in the production 
cross section will be reduced to 6-11%. The top mass uncertainty will be reduced to around 
1-2%. Clearly the agreement between standard model theory and experimental results is not 
close enough to include moderate shifts from the SM results. 

We have considered the complete one loop SUSY corrections to the process qq — > tt. 
These include SUSY-QCD and SUSY-QFD corrections not ignoring the box [both direct 
and crossed boxes]. Although box diagrams in general give a small contributions one must 
include them for completeness and exact numerical predictions. The purpose of this note is 
to concentrate on the complete one loop supersymmetric QCD corrections. Recently Li et. 
al. P] have reported the one loop SUSY-QCD corrections. We give a comparison between 
their and our work. In particular we find additional corrections [crossed box and gluon 
self-energy] which are not given by them. Several mistakes/misprints in their work are also 
noted, however their erratum [^ now corrects these. Importantly our numerical work does 
not agree with their original claim |^. However their numerical values given in the erratum 
|§] agrees now with the general trend that we give in this report. There still remains some 
disagreements as we report in section 4. We also give the energy dependence of the cross 
section and differential cross sections at the parton level. This allows us to compare our 
results with the one loop correction to the same subprocess in the standard model [0. One 
may use the parton level values to study the details of the SUSY corrections more directly. 

The one loop Electroweak corrections to the process qq — > tt have been considered 
by several groups [|10|iii the context of MSSM. Our results on these will be be presented in 
a subsequent paper. The complete SUSY corrections to the processes qq — > qq, qq — > qq, 
qq — > gg, and qg — > qg are being considered by [|TT|. 

The layout of this paper is as follows. In next section we give the one loop SQCD 
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radiative corrections to the process qq — > tt which arise from the gluon self-energy, the 
quark wave function renormahzation and the triangle diagrams. For completeness we also 
include the Born expression for the process qq — > tt. In Sec. 3, we write out the results for 
the corrections arising from the box diagrams [direct and crossed] due to the squarks and 
gluinos. Sec. 4 gives the numerical results. For our numerical work we use the Fortran code 



FF fl^for the evaluation of the scalar integrals [|l3l, and the MRSA parton distributions of 



Martin et. al. and finally the integrations are carried out by using BASES [|T5|. We have 
made several cross checks to make sure to eliminate any numerical errors. In the appendix 
we give the box contribution using the same momentum assignment of [§. We compare our 
results with |^ wherever required. 

II. TREE, AND THE ONE-LOOP CONTRIBUTIONS IN SQCD [EXCEPT FOR 

BOX ] TO THE PROCESS qq — > tt. 

At the parton level the processes responsible for the production of top[t] anti- 
top[f|in energetic pp collisions to order[a^ i.e. tree-level] are 

• The annihilation of quark-antiquark pair into top anti-top via a virtual gluon 
exchange 

q[pi]q[p2] — ^ t[p3]t[p4] 

• The fusion of gluon-pair into top anti-top via a virtual gluon or a virtual top- 
quark exchange 



g[Pl]g[P2] ^ t[p3]t[p4 
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Particle momenta have been shown in the parentheses. The schematic diagram for the first 
process is shown in Fig. 1, which is the reaction we choose to concentrate in this paper. 
To get a complete analysis one must include the second process as is done for the standard 
model [^], although it contributes only 10% at the Tevatron. We work with the Mandelstam 
variables s, t, and u defined as 



s= (P1 + P2)' = (P3 + P4)', (2) 

t = (Pl-P3)'=(P4-p2)', (3) 

u = (P2 - Pa)' = (P4 - Pi)' (4) 

The Mandlestam variables satisfy the relation s + 1 + u = 2m^ where we have taken the 
initial parton mass as zero. With our momentum assignments the leading order QCD matrix 
element of quark antiquark annihilation is given by 

iM^orn = ul(P3,S3)[-igsT^i7^]vUp4,S4)[-iV]v^(p2,S2)[-igsT^i7.]<(pi,Si) (5) 

s 

It is straightforward to obtain from the above equation the square of the Born matrix element 
averaged over initial spin and color degrees of freedom and summed over the final ones. We 
immediately obtain 

j:\MZrn\ = \Mo\' = ^F, (6) 

here and elsewhere in this paper we define 

Fi = 2sm? + (t-m,^)2 + (u-m,2)2 (7) 
The Born differential cross section is readily written as 



dt IGns' 
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Using the above equation to integrate over t, and noting that the integration hmits of t are 
in our case are given by the equation 

- ^§ - ^§/9t < t < - is + isA (9) 



[here = j^] immediately has for the expression for the Born cross section at 

the parton level 



„Born 
(7n = 



-Ms + 2ml] (10) 



3s2 ^9' 

We have intentionally written the above result in the form with the color factor separated 
out. The total amplitude squared upto one loop can be written as 

|Mo+l|' = |Mo|' + 2Re[MsewMS + MtnangleMj + MboxM^] (11) 

The total self-energy and wave-function renormalization [sew] contribution to the process 
qq — > tt can be written as 

M,ewMj = Uf^^^ul + M^qcdMJ (12) 

The gluon self energy diagram is shown in Fig. 2a. We renormalize the SUSY contribution 
to the QCD coupling as at zero momentum transfer. The gluon self-energy gets contribution 
from a gluino loop and s-quark loop. One may write 

Uf^^Ul = Mf mJ + Mf mJ (13) 
Mf = |Mop^(3) /'x(l-x)ln[l-x(l-x)4] (14) 

TT Jo 



MfMt = |Mor^(i) /' ^(l-2x)Mn[l-x(l-x)i (15) 

^ flavor ™qi 



The sum over two complex scalar fermions [16|, /i and /2, for each flavor is understood. For 



light quarks, the mass eigenstates are expected not to deviate significantly from the current 
eigenstates, fi and fn, due to chiral invariance. This may not be the case for the bottom 
and top squarks. To avoid plethora of indices we write all our results for the s-particle 1. 
However one must be careful that the total results do not always follow by replacing 1 by 2, 
for example in the box diagram we can have the mixed case i.e we may have squark [mass 
eigenstate 1] on one side and stop [mass eigenstate 2] on the other side of the box. We denote 
the squark mixing angle by 9 and that of stop by 9. The expressions for squark and stop 
mass eigenstates are 

qi = qL cos 6' + qR sin 6* (16) 
q2 = -qL sin 9 + cos 9 (17) 
ti = tL cos 6* + tR sin 6' (18) 

t2 = -tL sin 61 + tR cos 6' (19) 

The wavefunction renormalization [Fig. 2 a-d] contribution can be written as sum of 
two terms 

MsqcdMS = m"*mS + yrml (20) 

the [wt] contribution comes from the top wave-function renormalization and [wq] is from the 
parton [quark] wave function. 

yr'Ml = |Mo|2^[A2(Bi(m2,mg,m^) + 2mlB[) - B2(2mgmtB;)] (21) 

A? = 4 + hi (22) 



ai = — ^(cos^ — sin^) 
v2 



61 = — ^(cos 9 + sin 9) 
v2 



Aj = a? + b? 



and Bq to be used below is given by 

>2 _ ~2 1:2 



Oi = — ;= (cos 6^ — sin 9) 
v2 



61 = — ^ (cos ^ + sin 9) 
V2 

The total Triangle contribution to the process qq — > tt can be written as 

MtriangleMj = MtiMJ + MT2MJ3 + MtsMI, + Mt4MJ 

where [Fig. 3a] 



9s2 



MtiMS = ^VfTi 



Here i=l,2 



= ^['^rngmt{Bf)[Co + Cu] - 2m\{A\)[C2x + Cn] - 2{A\)C^^\F^ 



= ^[-2m-,mt{B^)[Co + Cn] + 2m'Mt)[C2i + Cn]]F2 
F2 = §2 (30) 

We note that Li et. al. P] have written 2sm^ + + s(s — 2m^) as the coefficient of their 
F5 which is equal to 2s^. In the above Eqs.^ and ^ the arguments of the C integral are 
Cij(-P3,P5,mg,mt- ,mt-). 

The expression for Mt2Mq [Fig. 3b] is rather simple. We can simply obtain it from 
the above by setting nit = 0. As a double check we have also calculated it directly and the 
result reads 

MT2Mt = M^Fr (31) 
ys J 

Here only i=l case is nonzero 

FP = ^[-2{Al)C,,]F, (32) 

In the above Eq.|32| the arguments of the C integral are Cij{—pi,p5,mg,mqj^,mq-^). 
We now give the expression for Fig. 3c 

^T3Ml = ^Y.Fr (33) 
ys . 

Here i=l,2 

= ^m-^^9MB?)[Co + Cn] + {A't)[[n - 2]C2, - s{C23 - C22) 

2,4:71 

- ml{Co + C21 + 2Cn) - mlCo]]F, (34) 



Fr = ^[9][2m,mtmCu + 2(A')K(Cn + C^^W^ (35) 
In the above Eqs.^ and |3^ the argument of the C integral are Cjj(— ^3, ps, ruf^ , nig, nig). 
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The triangle diagram for the qqg vertex is calculated directly and also as double check 
got from T3, by first replacing by ruq and then setting the latter equal to zero. One 
obtains 

MT4Mt = M^Fr (36) 

Here i=l 

= ^mm^ - 2]C24 - s{C23 - C,2) - mp,]]F, (37) 

In the above, Eq.|3^, the arguments of the C integral are Cij(— pi, ps, mq-^ , mg, mg). 

One can see from the above contributions of self energy ,wave function renormaliza- 
tion and triangles that they all factor into something times tree level amplitude except for 
contributions from triangle diagrams, Eqs.|3y and These arise since the top mass cannot 
be ignored! From the arguments of the above loop integrals we see immediately that they 
do not depend on the t-channel variable. From these simple observations one can see that 
the integration over t-channel variable for the above contributions is straightforward. This 
is not the case for the box diagrams since the box loop integrals depend explicitly on the t 
and u channel variables. 

III. CONTRIBUTION FROM THE BOX DIAGRAMS 

The total box contribution to the process qq — > tt can be written as 

MboxMS = + M^Z^l (38) 

where 

M?oM = 4^EfP^ (39) 

i 
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Here i=0,ll, 12, 13,23,24,25,26 and 27. By using the notation 

A+A+ = [a? + b?][a? + b?] + 4aibiaibi 

AsxAsx = [a? + hl][sil + hj] - 4aibiaibi 

A+A5X- [a?-b?][a? + b?] 
we find for the direct-box diagram [Fig. 4a] contribution: 

F^^^ = ^[m|(2A.A.)[2^m? + 2{t - m'.nD, 

TT 

F/^s ^ ^[^-^^(2A+A5.)[-252]]i^,, 

TT 

F,^^ = ^Km,(2A+A5,)[2s(S - 2m?) - 4(t - m?)^]]!?,^ 

TT 

F^^s = ^[-m5m*(2A+A5,)[2s(5 - 2m?) - 4(£ - m?)2]]A3 

TT 

F,^^ = ^K(2A+A+)[25m? + 2(t - m?)^ - 252]]D22 

TT 

^^2^ = ^[m?(2^+A+)[25m? + 2(t - m?)^ - 2s']]D,s 

TT 

F,^^ = ^K(2A+A+)[25^]]i^24 

TT 

= ^[(2^+^+)[25(^ - m?)2]]D,5 

TT 

= ^[(2A+A+)[2s^m? - 2m?(25m? + 2{i rnl f ) 2s{u - m?)^]]^: 

TT 
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TT 

In the above, the arguments of the D-functions are Di = Di[— pi, p3, p4, niq^, mg 
For the contribution of the crossed-box diagram [Fig. 4b] , we find 

M^o^Mt = [-^]^^FpB 

with 

A+A5x = -[a?+b?][&?-b?] 
F^^B ^ ^[ml{2A,,A,,)[2sm', + 2{u - m'.nDo 

TT 

rv 

FgB = ^[m~,mt{2AtA,,)[-2s']]Dn 

TT 

Fg"" = ^[mgmt(2AtA,,)[2s{s - 2m^) - A{u - m^,f]]Du 

TT 

pCB ^ '^[-mgmt{2AtA^,)[2s{s - 2ml) - 4(« - mlf]]D^^ 

TT 

Fg'' = ^[ml{2AtAt)[2sml + 2{u - mlf - 2f\\D,, 

TT 

FgB = ?^[ml{2AtAt)[2sml + 2{u - mlf - 2s^]]D,, 

TT 

Fg^ ^^[ml{2AtAt)[2s'']]D^ 

TT 

Fg,^^'^[{2AtAt)[2s{i-minD,, 

TT 

Ffe^ = —[{2AtAt)[2sV^ - 2ml{2sml + 2(a - m^,f) 2s{i - mlf]]D2. 

TT 
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FOB ^ ^[_2(2^+A+)[2Sm? + 2(t - m?)^]]D27 (64) 

In the above, the arguments of the D-functions are Di = Di[— pi, p4, pa, niq-^, mg, m^- , nig]. 

We are now in a position to write the expression for top pair production in proton anti- 
proton colhsion by weighing our expressions for differential cross section and cross section 
of the subprocess qq — > ti by the parton distribution functions and integrating over the 
parton variables i.e., 

da(pp ^tt) = Y. f dxi dx2[Di/p(xi, Q2) Dj/p(x2, Q')]da(ij ^ tt) (65) 

Here da represents the subprocess cross section at cm. energy square of s = X1X2S, where 
y/s is the cm. energy of the pp system. In our numerical calculation, we adopt the MRSA 



parametrization [|1^ for effective parton distribution evaluated at Q = . In order to show 



the SUSY radiative corrections clearly, we use a constant ag = 0.123. 



IV. NUMERICAL RESULTS AND CONCLUSIONS 

To facilitate comparison with work of ^ we give our numerical results for the same 
parameter values as in A more detailed numerical work will be given elsewhere. We thus 
take mt = 170 GeV, and assume no mixing between the squarks. The mass splitting between 
the squarks of different flavors is also ignored [§]. The common squark mass is denoted by 
mq. 

We first consider percentage one loop corrections at the hadronic cross section as a 
function of the squark mass. Taking the gluino mass to be 150 GeV we find the percentage 
corrections changes from —5.3% [mq = lOOGeV] to 1% [mq = 600GeV], see the solid curve 
of Fig. 5. This does not agree with the original claim of Li et al. |^., where they find 
for gluino mass of 150 GeV, 23% [mq = 100 GeV] and 5% [mq = 420 GeV]. However, the 
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corrected version |§[[see Erratum] values of —6% [niq = 100 GeVjand 4% [niq = 600 GeV] for 
gluino mass of 150 GeV are in more closer agreement with our values. The remaining small 
discrepancy may probably be explained since they p] have not included the gluon self-energy 
and the crossed box. For gluino mass of 200 GeV we find, see Fig. 5b, that the corrections 
change rapidly from 22% to —0.5% as squark mass changes from 100 GeV to 600 GeV. Here 
again there is no agreement with the original claim of ^ where they had reported a variation 
of 6.5% to 0% for a gluino mass of 200 GeV. The revised [erratum values of 31% [squark 
mass of 100 GeV] and 6% [squark mass 600 GeV] for gluino mass of 200 GeV, are still 
somewhat different from our values. We find for the gluino mass of 200 GeV, the relative 
corrections of 22%, 9%, 6%, 3%, 1%, -0.5% for mq = 100, 200, 300, 400, 500, 600 GeV 
respectively, which should be compared with 31%, 18%, 11%, 9%, 7%, 6% of the revised 
values of ||^. 

A comment is in order. It can be noticed from our Fig. 5 that the corrections change 
sign as the gluino mass is changed from 150 GeV to 200 GeV. As we have assumed a top 
mass of 170 GeV the threshold for top pair production is crossed in this region and hence 
the sign change and rapid change in magnitude of relative corrections occur. 

Next let us consider percentage one loop corrections at the hadronic differential cross 
section as a function of the squark mass. As an example, we show the case at the tt center- 
of-mass scattering angle 6cm = 10°. Taking the gluino mass to be 150 GeV we find the 
percentage corrections changes from —7.5% [mq = lOOGeV] to 2.5% [mq = 600GeV], see the 
solid curve in Fig. 6. For gluino mass of 200 GeV we find, see the dashed curve in Fig. 6, that 
the corrections change rapidly from 22% to 0% as squark mass changes from 100 GeV to 600 
GeV. As remarked before it is only the box loop correction which depends on the t-channel 
variable or on 6cm- If the box corrections are not larger than the other contributions, one 
would naively expect the percentage differential cross section to show only a weak dependence 
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on 9cm- This is indeed the case as can be seen by comparing Figs. 5 and 6. 

It is useful to give the percentage correction at the parton level since among other 
things they facilitate a comparison with correction found in the context of standard model 
0]. Moreover the corrections at the parton level can provide more direct and detailed tests 
of SUSY corrections. We first show in Fig. 7 the percentage corrections to the qq — > tt 
cross section as functions of the subprocess cm. energy At the parton level the total 
percentage cross section varies between 20% and —9% as center of mass energy is varied 
between 350 GeV to 1.8 TeV. When gluino mass is 150 GeV (the solid line), the correction 
is positive very near to the threshold (v^ < 430 GeV) but it gets rapidly negative down 
to -9% at around ~ 800 GeV. The resulting cancellation explains the smallness of 
the correction at the hadronic level; see the solid line at rriq = 200 GeV in Fig. 5. For 
rrig = 200 GeV (the dashed line), the correction grows from about 7% near the threshold to 
20% at = 400 GeV where the gluino-pair threshold opens. This large positive correction 
slighly above the ti threshold explains the large positive correction to the hadronic cross 
section in Figs. 5 and 6. We may compare these results to the standard model who 
report on the one loop virtual [electroweak] relative corrections to parton qq — > tt cross 
section among other things. As they take the top mass of 100 GeV and 250 GeV we can't 
compare our results directly with theirs. However we can extrapolate from their Figs. 9 and 
10 that for a top-quark mass of 170 GeV one would obtain corrections between 350 GeV and 
1.8 TeV of around 10% and -15%. 

In order to clarify structure of the SUSY radiative corrections, we show in Figs. 8a 
and 8b contributions from the gluon self-energy correction (dotted line), the sum of the 
triangle and quark wave-function corrections (short-dash), the direct-box (long-dash) and 
the crossed-box (dash-dotted) contributions separately. Fig. 8a is for rrig = 150 GeV and 
Fig. 8b is for nig = 200 GeV, while all the squark masses are set to 200 GeV in both 
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cases. While the gluon self-energy correction grows with energy, the sum of the triangle and 
quark wave-function corrections more than compensate for the effects. Direct box diagram 
contribution is positive near the threshold and turns to negative at higher energies. The 
crossed box diagram partially cancells the direct box contribution. Box contributions are 
generally smaller than the other corrections as expected from the similarity of the corrections 
to the total and differential cross sections. The gluino-pair threshold effects are evident in 
all the curves in Fig. 8b. 

We finally consider the percentage corrections at the parton level taking -\/§ = 600 GeV 
and letting the squark mass vary between 100 GeV and 600 GeV. For the gluino mass of 150 
GeV we find that the percentage corrections of differential cross section vary between —20% 
and 6.5%, [Fig. 9]. For the gluino mass of 200 GeV the corrections vary between —7% and 
4%, [Fig. 9]. 

Note added:After the calculation was completed and the present paper was being 
written up: the following works came to our attention: 



1: J. Kim et al. examine both the SUSY Electroweak and SUSY QCD like correc- 
tion. However they do not include box diagrams and claim that the box contributions are 



small citing J. Ellis and D. Ross, [|TT| and P. Kraus and F. Wilczek, [|T^ works as evidence. 
These authors state that their results for SUSY QCD agree with Ref. |P| while those of SUSY 
Electroweak disagree with J. Yang and C.S. Li [|T0|. 



2: J. Ellis and D. Ross [^ work at the parton level considering the processes qq — > qq 
qq — > qq, qq — > gg, and qg — > qg. However they do not consider tt cross section as it 
requires separate treatment. 



3: P. Krauss and F. Wilczek [0 also studied the SUSY corrections to the quark gluon 
scattering processes in the limit of large SUSY particle masses. 
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FIGURES 

Fig. 1 Tree-Level Diagram for the process qq — > ti. 

Fig. 2a Schematic diagram for the Gluon Self-Energy due to SQCD particles. 

Fig. 2b Schematic diagram for the Quark WFNR due to the SQCD particles. 

Fig. 2c Schematic diagram for the Anti-Quark WFNR due to the SQCD particles. 

Fig. 2d Schematic diagram for the Top WFNR due to the SQCD particles. 

Fig. 2e Schematic diagram for the Anti-Top WFNR due to the SQCD particles. 

Fig. 3a Triangle contribution from two stops and one gluino to the ttg vertex. 

Fig. 3b Triangle contribution from two squarks and one gluino to the qqg vertex. 

Fig. 3c Triangle contribution from one stop and two gluinos to the ttg vertex. 

Fig. 3d Triangle contribution from one squark and two gluinos to the qqg vertex. 

Fig. 4a Direct Box contribution from one stop, one squark and two gluinos. 

Fig. 4b Crossed Box contribution from one stop, one squark and two gluinos. 
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Fig. 5 One-loop percentage relative corrections to the tt production cross section in pp colli- 
sions at = 1.8 TeV as a function of squark mass for a gluino mass of 150 GeV (solid 
line) and 200 GeV (dashed line). The top quark mass is set to 170 GeV. 

Fig. 6 One-loop percentage relative corrections to the tt production differential cross section 
at 6'cm = 10° in pp coUisions at = 1.8 TeV. The results are shown for mt = 170 GeV 
as functions of the squark mass for a gluino mass of 150 GeV (solid hne) and 200 GeV 
(dashed line). 

Fig. 7 One-loop percentage relative correction to the parton-level qq — > tt cross section for 
mt = 170 GeV as functions of the qq center of mass energy, -\/§. The results are shown 
for mq = 200 GeV and mg = 150 GeV (sohd line) or m^ = 200 GeV (dashed line). 

Fig. 8a Individual contributions of the one-loop percentage relative correction to the qq — > tt 
cross section are given as functions of the qq cm. energy v^. Gluon self-energy 
correction (dotted line), the sum of triangle and quark wave- function corrections (short- 
dashed hne), the direct-box diagram (long-dashed hne) and the crossed-box diagram 
contributions (dash-dotted line) are shown separately. Also shown by the solid line is 
the total percentage relative correction, as given in Fig. 7. The mass paramters are: 
mt = 170 GeV, ruq = 200 GeV, and m~g = 150 GeV. 

Fig. 8b Same as Fig. 8a except rrig — 200 GeV. 

Fig. 9 One-loop percentage relative corrections to the parton-level qq — > tt differential cross 
section at ^cm = 10° and — 600 GeV. The results are shown for mt = 170 GeV as 
functions of the squark mass for a gluino mass of 150 GeV (solid hne) and 200 GeV 
(dashed line). 
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APPENDIX: CONTRIBUTION FROM THE BOX DIAGRAMS 



As already mentioned we give in this appendix the box results in the notation of 
Li et al. for the purposes of exact comparison. The total box contribution to the process 
qq — > tt can be written as 

MboxMj = MZmI + M'^^M (Al) 

where 

m?oM = ^EfP^ (A2) 

Here 1=0,11,12,13,23,24,25,26 and 27. 

A^AjI" = [ai + h\] [ai + b^] + 4aibiaibi (A3) 

As.As, = [a^ + b^] [&2 + b^] - 4aibi&ibi (A4) 

A+A5.= [a?-b?][&? + b?] (A5) 
We note that the following relations hold between our A's and cr's of |^ 

2A5..A5. = 4 + (A6) 



Fo^^ = ^[m? (2^5,^5x0 [2sm2 + 2(t - m?)^] + mgmt{2AtA,,)[2s^]]Do (A7) 

77 



F^^^ = ^[mgmt{2AtA,,)[2s']]Du (A8) 
vr 

21 



^DB ^ l!£[_^2^2A+A+)[2sm? + 2{u - m?)^] - m-,mti2AtA,,)[2s']]Di2 (A9) 

TC 

TT 

-mgmt{2AtA,^)[2sis-2ml)-A{i-m^,y]]D,, (AlO) 
F2T = -H{2AtAt)[2sm', + 2{i - m\f - 2s^\\D2, (All) 

TT 

F^^i^ = ^[-(2A+A+)[25^m? + 2S(^ - mlfWD^, (A12) 

TT 

F,^^ = '^\{2AtAX){2s{u - m\nD,, (A13) 

TT 

^^2^ = -K^(2^^^^)[2S^]]I^26 (A14) 

TT 

F^^^ = '^\-2{2A\A\)\2sm\ + 2(n - m,^)^]]^^^ (A15) 



Here the arguments of the D-functions are Dj = Di[— pi, — P2, P4, nig, niq. , nig, mfj. The above 
results for the direct box agree with |]^ after taking account of their erratum. 

Although Ref. @] does not give the contributions of the crossed-box diagram, we give 
our result in their notation. 

mM = [-^]^EfP (A16) 

Here 1=0,11,12,13,23,24,25,26 and 27. We note that the crossed box color factor is smaller by 
a factor of | compared to the direct box, and it contributes destructively with the direct-box 
contribution; see Eq. |AT| . 

A+A^=-[a? + b?][a?-b?] (A17) 
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pCB = ^[ml{2A,,A,,)[2sml + 2{u - m?)^] + mgmt{2At A,,)[2s^]]Do (A18) 

TT 

F^B ^ ^[rngmt{2A[A^,)[2s']]D^, (A19) 

TT 

pCB ^ ^[_^2^2A+A+)[2sm2 + 2{i-mlf] - mgmt{2At A,,)[2s^]]Di2 (A20) 

TT 

Fg^ = ^[mU2AtAt)[2{i-m',r] 

TT 

- m^mt(2]4p^)[2s(s - 2771^) - 4(ti - m^j^jjDia (A21) 
Fg^ = ^[m?(2A+A+)[25m,2 + 2{u - m\f - 2s^\\D^, (A22) 

TT 

i^CB ^ ^[_(2A+A^-)[25 V + 2^(i - m,^)2]]D24 (A23) 

TT 

F,^^ = ^[(2A+A+)[2s(^ - mlfWD^, (A24) 

TT 

i^2^^ = -\m\{2AlAl)\2n^^, (A25) 

TT 

F,^/ = '^{-2{2 At At){2sm\ + 2(£ - m,^)2]]i?27 (A26) 



Here the arguments of the D-functions are Di = Di[— pi, — p2, Ps, nig, niq. , nig, mfj. 
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Fig.l: Tree-level Diagram 




Fig. 2a: Gluon Self Energy due to Particles 



Fig. 2b: Quark Self Energy due to SQCD Particles 
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Fig. 2c Anti-Quark Self Energy due to SQCD Particles 




Fig.2d Top Self Energy due to SQCD Particles 




Fig.2e Anti-Top Self Energy due to SQCD Particles 
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Fig. 3a: Triangle contribution from two stops and one gluino 



Fig. 3b: Triangle contribution from two squarks and one gluino 




Fig. 3c: Triangle contribution from two gluinos and one stops 
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Fig. 3d: Triangle contribution from two gluinos and one squark 



Fig.4a: Direct Box contribution from SQCD Particles 




Fig.4b:Crossed Box contribution from SQCD Particles 



27 



00 



30 



20 



< 







■10 



100 



Top quark mass = 1 70 GeV 

Gluino mass = 150 GeV 

Gluino mass = 200 GeV 



300 400 500 

Squark mass [GeV] 



30 



20 



10 








■10 



500 



700 



Top quark mass = 1 70 GeV 

Squark mass = 200 GeV 

Gluino mass = 150 GeV 
Gluino mass= 200 GeV 



900 1100 1300 
Vs in GeV 



1500 1700 



Fig. 7 



CO 



o 

CO 

c 
o 

o 

CD 

o 
O 



30 
25 
20 
15 
10 
5 

-5 
-10 



Crossed Box 
Direct Box 
Gluon Self Energy 
Traingle+Quark w.f. 
Total 




Top quark mass = 1 70 GeV 
Gluino mass = 150 GeV 
Squark mass = 200 GeV 



■15 



-20 
-25 



500 



1000 
Vs in GeV 



1500 



Fig. 8a 



CO 



1^^ 

A 
ID- 

o; 
<D 

O 

w 
c 
o 

'■4— ' 

o 

o 
O 



30 
25 
20 
15 



5 

-5 
-10 
-15 
-20 
-25 



7 (1 Gluon Self Energy 
I Triangle+ Quark w.f. 

- , ^ Crossed Box 

1 \ \ Direct Box 


II 

Tnn niiark ma^^ 
Gluino mass = 
Squark mass 


1 1 

- 1 70 GeV 
200 GeV 
= 200 GeV 


- , \ \ Total 


















^^^^^^^ ^ 


— 




- 


" - ^ 

II 


III" 


500 


1000 
Vs in GeV 

Fig. 8b 


1500 




-30 

100 200 300 400 500 600 

Squark mass [GeV] 
Fig. 9 



